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Abstract 


The acoustics of an advanced single rotation 
SR-3 propeller at cruise conditions are studied 
employing a time-domain approach. The study evalu- 
ates the acoustic significance of the differences 
in blade pressures computed using nonreflecting 
rather than hard wall boundary conditions in the 
three-dimensional Euler code solution. The direc- 
tivities of the harmonics of the blade passing fre- 
quency tone and the effects of chordwise loading on 
tone directivity are examined. The results show 
that the maximum difference in the computed sound 
pressure levels due to the use of blade pressure 
distributions obtained with the nonreflecting 
rather than the hard wall boundary conditions is 
about 1.5 dB. The blade passing frequency tone 
directivity obtained in the present study shows 
good agreement with Jetstar flight data. 

Introduction 


Advanced high speed propellers incorporate 
many features which are quite different from those 
of conventional low speed propellers. Advanced 
propellers employ eight or ten highly loaded, low 
aspect ratio, thin and highly swept blades. These 
high speed propellers with transonic helical tip 
Mach numbers produce significantly more noise than 
low speed propellers. Therefore, special attention 
must be paid to reducing high speed propeller noise. 
To accomplish this, an understanding of the noise 
generation and radiation mechanisms and noise con- 
trol techniques is essential. Analytical models, 
numerical approaches and experimental methods are 
being developed to understand and predict high- 
speed propeller noise levels. I -4 


Propeller noise consists of two dominant com- 
ponents: (1) Thickness noise and {2) loading 

noise. The thickness noise component is governed 
by the geometry of the blade and can be computed 
reasonably accurately. The loading noise component 
depends on the blade pressure distribution, which 
has often been estimated by approximate methods for 
input to the acoustic computations. With the 
availability of super-computers the blade pressure 
distribution can be obtained more accurately by 
solving the full three-dimensional Euler equations 
governing the flow field of a high speed propeller. 

The blade pressure distributions obtained from 
the Euler analysis of the propeller flow field can 
be employed to compute the sound pressure levels of 
advanced propeller tones by a procedure such as 
that used by Clark and Scott. 5 They modified the 


Denton code 5 developed for the three-dimensional 
analysis of the flow in turbomachines to compute 
the propeller flow field. Their modification con- 
sisted of moving the casing away from the blade tip 
to a distance of twice the blade radius from the 
axis of the propeller. In other words, the flow 
field of the propeller, which is unconfined, was 
modeled as though the propeller was inside a much 
larger solid-walled duct. It has been arqued that 
the hard wall boundary condition could signifi- 
cantly affect the computed blade pressure distribu- 
tions, and the noise levels computed from these 
pressure distributions could be significantly in 
error. 

The present investigation was undertaken to 
establish the influence of the far field boundary 
condition on the blade pressure distribution and on 
the noise levels computed from them. The hard wall 
boundary condition is replaced by a nonref lectinq 
boundary condition. The characteristic variable 
corresponding to the incoming acoustic information 
is imposed on the boundary. The blade pressure 
distributions obtained with this boundary condition 
are compared with those obtained with the hard wall 
boundary condition. The acoustic significance of 
the change in blade pressure distributions due to 
the implementation of the nonreflecting boundary 
conditions is evaluated. This is done by comput- 
ing blade passing tone directivities using the 
blade pressure distributions obtained with the two 
boundary conditions as input to Farassat's computer 
program 7 for sound pressure levels. The computed 
directivities are compared with flight data 
obtained on the SR-3 model propeller installed on 
the Jetstar test bed aircraft as shown in Fiq. 1. 
The computed directivities of the higher harmonics 
of the blade passing frequency and the effect of 
chordwise loading on directivity are also examined. 

Computation of Blade Pressure Distributions 


In the present study we consider the acoustics 
of the SR-3 high-speed propeller at design condi- 
tions. The blade pressure distributions are 
obtained by solving the three-dimensional Euler 
equations employing Denton's time marching method. 
As mentioned above, the hard wall boundary condi- 
tion of Ref. 5 for the outer boundary of the compu- 
tational domain representing the propeller flow 
field is replaced by a nonreflecting boundary con- 
dition. The specification of the nonreflectinq 
boundary condition reduces the influence of the 
location of the far field boundary on the solution. 
Further details of the nonreflecting boundary con- 
dition employed, the aerodynamic solutions 
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obtained, and the comparison of the solutions with 
data and other predictions are reported in Ref. 8. 
The blade pressure distributions obtained for the 
SR-3 propeller blade at design conditions (advance 
ratio of 3.06 and a flight Mach number of 0.8) are 
shown in Fig. 2. Pressure distributions on the 
blade are shown at 21, 53, and 84 percent span. 
First, it is seen that significant differences 
exist between the pressure distributions obtained 
with the hard wall and the nonreflecting boundary 
conditions. The higher pressures with the hard 
wall boundary condition are attributed to the pres- 
sure rise due to the confining wall. It is also 
shown in Ref. 8 that the results obtained with the 
nonreflecting boundary condition are in better 
agreement with the experimental data than those 
obtained with a hard wall boundary condition. The 
computed total power coefficient, elemental power 
coefficient, elemental thrust coefficient and the 
flow swirl angle downstream of the blade were com- 
pared with the experimental data. 


where c and p Q are the speed of sound and 

density in the undisturbed medium respectively, 
and V n is the local normal velocity on the 
blade surface, which is described by the equation 
f(x,t) - 0. The local force (per unit area) on the 
fluid by the blade surface is 1 -j , and e(f) is 
the Dirac delta function. The following two forms 
of the solution of Eq. (1) are used to calculate 
propeller noise in Farassat’s code: 


4irp’(x,t) 


1 a_ 
C 3 1 



ds 

Iret 



(subsonic) 
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Figure 3 shows chordwise loading distributions 
in terms of normalized differential pressure AC p 

at seven spanwise locations, namely 16, 26, 37, 53, 
63, 79, and 95 percent span. With the nonreflecting 
boundary condition, the pressure difference near 
the leading edge is reduced and the peak near the 
trailing edge is sharpened compared to that for the 
hard wall boundary condition. This chordwise load- 
ing distribution influences the directivity of the 
blade passing frequency tone as discussed in the 
Results section. The peak in the loading distribu- 
tion near the trailing edge is due to the trailing 
edge shock system observed in the aerodynamic anal- 
ysis and experiments.* It is instructive to note 
here that the approximate methods used to calculate 
the loading distribution, such as that of Hanson 9 
(see Fig. 14 of Ref. 9), show little evidence of a 
trailing edge shock system. The radial distribu- 
tions of the loading for the two boundary condi- 
tions are shown in Fig. 4. The hard wall boundary 
condition results in a higher loading than the non- 
reflecting boundary condition. The power coeffi- 
cient calculated with the hard wall boundary 
condition is about 11 percent higher than the 
measured value at design conditions, while that 
obtained with the nonref lecting boundary condition 
is about 3 percent lower than the measured value. 

Computation of Advanced Propeller Noise 

The blade pressure distributions obtained 
above are used as input to the acoustic computation 
employing Farassat's code.?’ 10 This computer 
code for advanced propeller noise prediction uses a 
time domain formulation. It is valid for both near 
and far field noise calculations and handles 
observers fixed to the ground frame or fixed to the 
aircraft frame. His formulation uses two forms of 
the solution of the Ffowcs-Williams and Hawking 
equation with thickness and loading source terms 
only. 


, C p„CV„ + 1 
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where ds is an element of surface area of the 
blade f - 0 and dr is an element of length of 
the curve of intersection for the surface f = 0 and 
the radiation vector r = x - y. Suffix r refers 

to the radiation direction. M*. = v . r/c, where 
v is the local velocity on the blade surface and 

r is the unit vector in the radiation direction, 
e is the angle between the surface normal and 
the radiation vector, g = t - t + r/c =0 is 
the equation of the collapsing sphere centered at 
the observer position x, t is the observer time, 
and t is the source time. Equation (2) is used for 
1^ < 0.98 and Eq. (3) is used for M r > 0.98. The 
equations are solved numerically using finite 
difference techniques. 


The upper and lower surfaces of the blade are 
divided into panels. The computed pressure distri- 
bution (from the aerodynamic code) on these panels 
is used to calculate the acoustic pressure signa- 
ture at any observer position. The contributions 
of all panels are added together to obtain the 
acoustic pressure p'(x,t). The pressure signature 
for only one blade is computed. The signature for 
several blades is calculated by shifting the signa- 
ture for one blade in time and summing the pres- 
sures for each observer time within a period. The 
acoustic pressure signature is then Fourier ana- 
lyzed to obtain the spectrum of noise. 

Discussion of Results 


The governing equation of the acoustic pres- 
sure p'(x,t) of a propeller? is 


1 3 p ' 

2 2 
c^ 3 1 


? 

rp' 



( 1 ) 


The computations have been carried out for the 
eight bladed SR-3 model propeller configuration for 
the design conditions - flight Mach number of 0.8, 
flight altitude of 35 000 ft and an advance ratio 
of 3.06. The detailed pressure distributions on 
the blade surfaces obtained from the aerodynamic 
calculations are used as input to the acoustic com- 
putations. The results of the computations and 
comparisons with the Jetstar flight data 11 are 


presented and discussed. The microphone positions 
in the flight test are shown in Fig. 5 to aid 
discussion of the results. 

Boundary Condition Effects 

First we examine the effect of the far field 
boundary condition employed in the aerodynamic 
calculations on the acoustic results. The blade 
pressure distributions obtained with the hard wall 
boundary condition and the nonreflecting boundary 
condition at the outer boundary of the computa- 
tional domain for aerodynamic calculations are used 
to calculate the sound pressure levels of the SR-3 
propeller. Figure 6 shows the directivity of the 
loading noise component of the blade passing fre- 
quency (BPF) tone level for the two cases. In this 
figure and in the subsequent ones for directivity, 
the "predicted" curves are lines drawn through all 
the predicted points. Sound pressure levels com- 
puted with the pressure distributions obtained 
with the nonreflecting boundary condition show an 
increase in SPL ahead of the peak location and a 
small decrease behind it, compared to those 
obtained with the hard wall boundary condition. 

The maximum difference in the SPL occurs in the 
plane of the propeller and is about 2 dB. It is 
also seen that a local minimum exists in the plane 
of the propeller in both cases. 

The overall sound pressure levels (loading 
plus thickness) are shown in Fig. 7. Also shown 
are the Jetstar flight test data 11 for comparison. 
The flight data were corrected for the effects of 
refraction in the fuselage boundary layer and scat- 
tering at the surface of the fuselage or boom based 
on Ref. 12 (see also Ref. 5). The trends of the 
two curves are nearly the same as in Fig. 6: in 
the nonreflecting boundary condition case the pre- 
dicted SPL is higher ahead of the peak location and 
slightly lower behind it, compared to the hard wall 
boundary condition case. The maximum difference 
between the two cases is at the plane of the pro- 
peller and is about 1.5 dB. It is interesting to 
note that although the far field hard wall boundary 
condition is found to result in significant over- 
predictions in the aerodynamic calculations, 8 it 
affects the predicted SPL only marginally. 

The predicted sound pressure levels agree 
fairly well with data, with the predictions using 
the nonreflecting boundary condition giving a 
slightly better agreement. The reason for the 
rapid fall -off of the fuselage microphone data in 
the aft quadrant is not known. One point to note 
here is the following: although the loading 

obtained with the nonreflecting boundary condition 
is less than that with the hard wall boundary con- 
dition (Fig. 4), the SPL predicted ahead of the 
propeller is higher in the nonreflecting case. 

This is due to the differences in chordwise loading 
distribution between the two cases and will be 
discussed in more detail later in this paper. 

Figure 8 shows the directivity of the second 
harmonic of the BPF. The predicted loading noise 
component and the total noise level for the two 
boundary conditions are shown. The trends of the 
two curves are similar to those observed in 
Figs. 6 and 7. However, the differences between 
the SPL levels for the two boundary conditions are 
small. The agreement of the predicted total noise 
level with the data is good. 


The spectrum of the loading component of the 
propeller noise is shown in Fig. 9, for the two 
boundary conditions. Noticeable differences in SPL 
levels of the two cases occur only for the first 
and second harmonics. Figure 10 shows the spectrum 
of overall noise in the plane of the propeller. 
Again, the boundary conditions affect the predicted 
SPL levels of only the first two harmonics. 

Further discussions will be based on the 
results obtained using the blade pressure distri- 
butions obtained with the nonreflecting boundary 
condition. 

Chordwise Loading and Directivity 

Efforts to understand the directivity obtained 
with the pressure distribution for a nonreflecting 
boundary condition (Figs. 6 and 7) led to an exami- 
nation of the effect of chordwise loading on direc- 
tivity. The fact that accurate noise predictions 
for the propeller require accurate estimates of the 
chordwise loading distribution has been stressed by 
Hanson-*- 3 and Farassat.? Hanson observes that 
(1) "peaked" chordwise loading distributions should 
be avoided to take advantage of the chordwise non- 
compactness, and (2) shifting the loading to the 
trailing edge causes a phase lag which would be 
used to promote phase cancellation between signals 
from different radii of the blade. However, the 
effect of these two factors on directivity is not 
clear. 

To investigate the effects of chordwise load- 
ing distribution on the tone directivity, the load- 
ing was increased near the leading edge and 
decreased near the trailing edge along the entire 
span of the blade. Figures 11(a) and (b) show the 
modified chordwise (at 37 percent span, as typical) 
and spanwise loading distributions. This modified 
loading resulted in a 9 percent increase in blade 
power. The directivities obtained with these load- 
ings are compared in Fig. 11(c) at the microphone 
positions of flight test. The sound pressure 
levels in these comparisons have been scaled to the 
design power. The modified loading distribution 
results in lower tone levels forward of the pro- 
peller plane. The directivity effect is similar 
to the results with the hard wall boundary condi- 
tion in Figs. 6 to 8, when a correction is applied 
for the difference in predicted total mechanical 
power levels. 

Acoustic Pressure Signatures and Spectra 

The computed acoustic pressure signature in 
the plane of the propeller (corresponding to the 
boom microphone 3) is shown in Fig. 12 along with 
the measured data. 34 ’ 38 The agreement of the 
predicted signature with data is good. 

Farassat 34 * 38 observes that the shape of the 
acoustic pressure signature depends on the chord- 
wise loading distribution. He finds that a tri- 
angular chordwise loading with the peak at the 
leading edge produces a wave form 38 which is in 
good agreement with data. 

The predicted and measured spectra correspond- 
ing to the boom microphone positions are shown in 
Fig. 13. The agreement of the predicted levels 
with the data is good for the first three harmonics. 
The reason for the discrepancies at higher 
harmonics are not known. The new formulation of 
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Farassat 14 * 16 also shows such discrepancies. 
Frequency domain calculations of higher harmonics 
of advanced propellers have not been reported. For 
low speed propellers, frequency domain calculations 
of Williams 1 ' and time domain calculations of 
Farassat 18 show good agreement with the data up to 
six harmonics. 

Directivities of the BPF Harmonics 

Here we examine the directivities of the har- 
monics of the BPF and also those of the individual 
loading and thickness components of the harmonics. 
Figure 14(a) shows the thickness , loading, and 
total noise directivities of the blade passing fre- 
quency tone level. The thickness noise shows a 
maximum just aft of the plane of rotation of the 
propeller. The loading noise component shows a 
local minimum in the plane of the propeller. The 
directivity of the total noise shows that the maxi- 
mum noise level occurs aft of the propeller at 
about 0.25 meters from the plane of rotation, which 
coincides with the location indicated by the flight 
data. The frequency domain calculations of the 
directivity of the total noise reported in 
Refs. 12 and 15 also show the maximum to occur aft 
of the plane of rotation of the propeller, but 
indicate a local minimum to occur in the plane of 
rotation. These calculations employ the near field 
frequency domain theory of Hanson. 1 ^ The total 
noise directivity calculations reported in Ref. 20, 
which uses the far field frequency domain theory of 
Hanson 13 *^ 1 does not show a local minimum in the 
plane of the propeller. The differences between 
the calculations of the time domain and frequency 
domain 1 ^* 18 approaches could be due to the differ- 
ences in input loading distributions. 

Figure 14(b) shows the directivities of the 
second harmonics and its components. The total 
noise directivity shows a fairly good agreement 
with data. Similar agreement of the total noise 
directivity with the data is also observed for the 
third harmonics (Fig. 14(c)). The predicted total 
noise of the fourth and fifth harmonics (Figs. 14(d) 
and (e)) are significantly different from the data. 
For these and higher harmonics the present numeri- 
cal method appears inadequate. Figure 15 shows the 
directi vities of the total noise of the first five 
harmonics of the BPF for comparison. As the har- 
monic number increases, the SPL at all angles tends 
to decrease. Amplitude irregularities around the 
plane of rotation at the fourth and fifth harmonics 
may be artifacts of the numerics. Detailed experi- 
mental directivities are unavailable. 

Conclusions 

The acoustic significance of the far field 
boundary conditions (hard wall and nonreflecting) 
employed in aerodynamic calculations has been 
studied. The present computations show that the 
maximum difference in the computed SPL due to the 
use of a hard wall boundary condition, in place of 
a nonreflecting boundary condition, is about 1.5 
dB. While the acoustic difference is small, the 
nonreflecting boundary condition is clearly more 
correct physically, and has been shown to give 
superior aerodynamic results. The location and 
level of the maximum of the BPF tone level pre- 
dicted by the present calculations are in agreement 
with flight data. The chordwise loading distribu- 


tion influences the relative levels of the SPL fore 
and aft of the propeller. It appears that the 
prediction of harmonics higher than the third needs 
refined approaches. 
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FIGURE 1.- SR-3 PROPELLER INSTALLED ON JETSTAR AIRCRAFT FOR FLIGHT TEST. 
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FIGURE 5.- NICROPHONE POSITIONS IN JETSTAR FLIGHT TEST. 
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FIGURE 6.- CALCULATED DIRECTIVITIES OF LOADING NOISE 
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FIGURE 7.- DIRECTIVITY OT TOTAL NOISE AT BPF (DESIGN- 
7348. 5 RPM; RUN 40 = 7314 RPM; RUN 42 = 7569 RPM). 
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FIGURE 8.- DIRECTIVITY OF 2 X BPF. 
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FIGURE 9.- CALCULATED LOADING SPECTRA AT 0.81 M SIDE- 
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FIGURE 10. COMPUTED TOTAL SPECTRA AT 0.81 M SIDELINE TO 
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FIGURE 11.- EFFECT OF CHORDWISE LOADING ON BPF DIRECTIVITY. 
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(B) EXPERIMENT. 

FIGURE 12.- WAVEFORMS IN THE PLANE OF THE PROPELLER 
(BOOM MICROPHONE 3). 





(A) NOISE SPECTRUM IN THE PLANE OF THE PROPELLER: BOOM 
MICROPHONE 3. 


FIGURE 13.- NOISE SPECTRA. 
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(A) DIRECTIVITIES OF BPF. 

FIGURE 14. - DIRECTIVITIES OF BPF AND ITS HARMONICS. 
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FIGURE 14.- CONTINUED. 
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